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PARTON DISTRIBUTIONS DIS06. 
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I discuss the current status of parton distributions. I outline the wide variety of 
different parton distributions available, and highlight which are either necessary or 
suitable for use at present. 

There are a large number of different parton distributions. If we consider 
the different types of particle which are partons we start with the quark 
model valence partons Q^) and dvix, Q^)- However, these carry only 

50% of the proton's momentum: there are also u{x, Q^) and d{x, Q^), which 
are not the same^; gluons g{x,Q^) which carry over 30% of the momen- 
tum; s{x,Q'^) and s{x,Q'^), with the possibihty^ that s{x,Q'^) ^ s{x,Q'^); 
c{x,Q'^) and &(a;, Q^), which are perturbatively generated (there could be 
intrinsic contributions'^ at large x); and at some level isospin violation, 
i.e. uP{x, Q^) 7^ d"-{x, Q"^), uP{x, Q'^) ^ d"-{x, Q^) - automatically the case 
with QED corrections'*, which also lead to a 7(0;, Q^) distribution. Over- 
all there are 6 — 16 different parton distributions, but some are very small 
and are often not needed. There is another way of counting, i.e. the dif- 
ferent sets of parton distributions from different prescriptions: LO, NLO, 
or NNLO in as', with resummation corrections or allowances for higher 
twist; using MS, DIS, or potentially other factorization schemes; fixed- 
flavour (FFNS), zero-mass variable-flavour (ZM-VFNS) and general-mass 
variable- flavour number scheme GM-VFNS, and even different versions of 
the last. This freedom, and choices in data sets used, cuts apphed, ways 
to treat errors, etc. lead to a staggering array - CTEQ4A1, CTEQ4A2^ ... 
CTEQ5HJ6, ... CTEQ6^ MRST98** ... MRSTOSc^, ... MRST04QED4, 
MRST04i°..., AlekhinOO", AlckhinOSi^, GRV98i3, Fermi02"..., ZEUSl^ 
ZEUS-ZJ^^, Hl*^, Botjc*^ and many more. Are all of these different sets 
really necessary? This is a complicated and controversial question. 

One restriction is very easy to impose - many (still used) partons are 
simply out of date. Unless there is a very good reason, one should not 
use pre-2000 parton distributions. The available data have improved a 
great deal since then, particularly the HERA structure functions*^' and 
Tevatron jets^". Also, some older partons have minor bugs. 
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Figure 1. Comparison of the NLO up distribution with the NNLO up distribution. 

Consider the order of the partons. LO requires as{M'^) ^ 0.130, and 
the from a global fit is far inferior to that at NLO and NNLO^^. The 
LO partons are qualitatively different from NLO and NNLO partons in 
some regions, in particular the gluon is much bigger at small x due to 
important corrections in the splitting functions. This can cause misleading 
conclusions on the evidence for saturation etc.. Such results derived from 
LO partons should be treated with care. The default has long been NLO, 
but the NNLO coefficient functions for structure functions have long been 
known^^ and the splitting functions are now complete^'^. These improve 
the quality of the fit slightly^^ and reduce as. A big change in the partons 
can occur when going from NLO NNLO, as seen in Fig. 1. 

To perform an absolutely correct NNLO fit we need both exact NNLO 
splitting functions and cross-sections. The NNLO Drell-Yan cross-sections 
have recently been calculated as a function of rapidity^^, leading to a de- 
crease in the sea quarks. The one remaining gap is the NNLO corrections to 
jet production in pp{p) collisions. However, the NLO corrections themselves 
are not large - at central rapidities they are < 10%, similar to the size of 
the correlated errors. There are also some NNLO estimates, i.e. the leading 
threshold corrections, which are expected to be a significant component of 
the total^^ (there are issues concerning the application within a given jet 
definition). These give a flat, small 3 — 4% correction, which is consistent 
with what we already see at NLO and is much smaller than the systematic 
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errors on the data. Hence, it seems perverse to leave the jet data out of a 
NNLO fit due to the lack of the full NNLO hard cross-section. 

For a full NNLO fit we also require a rigorous treatment of heavy quark 
thresholds, which is now available^^. Therefore, an essentially full NNLO 
determination of partons is possible. Of course, this is the best way to test 
our understanding of the partons in terms of QCD, but we only know a lim- 
ited number of cross-sections at NNLO. Processes with two strongly inter- 
acting particles are largely completed - DIS coefficient functions, pp{p) — > 
j*,W,Z (including the rapidity distribution), H,A'^,WH,ZH. For many 
other final states the NNLO cross-sections are not known and NLO is still 
more appropriate. Moreover, resummations may be important even beyond 
NNLO in some regions, as may higher twist. 

There is the issue of factorization schemes. In practice cross-sections 
are calculated in MS scheme, so we use MS parton distributions. However, 
DIS-scheme can be more useful for relating partons to real physical results, 
e.g. it is easier to fit the Tevatron jet data in the DIS scheme^^, or viewed 
differently the large high- a; gluon required at NLO and NNLO in MS scheme 
can be determined from the scheme dependence^". Schemes other than MS 
are valuable in this type of context, but at present rarely used. 

There are also partons corresponding to different prescriptions for heavy 
flavours. FFNS is intrinsically inferior to VFNS - it does not sum In /m\ 
terms in the perturbative expansion and at high scales this can lead to inac- 
curacies. Moreover, it is often necessary to have heavy flavour partons due 
to the lack of mass effects in the known cross-sections. Nevertheless, FFNS 
partons are also sometimes needed because some hard cross-sections are 
only calculated in this scheme^^'^^'^'^. However, in this case the treatment 
must be correct, and is often not so. The NLO (0(a|)) coefficient functions 
for heavy flavour in DIS are calculated^^ in a renormalization scheme where 
the coupling ag is fixed at 3 flavours. The partons have to be defined in 
the same way, otherwise there is double counting of q;| In {Q'^ /rnj^), terms 
which can lead to large errors^^. Also, there are no FFNS coefficient func- 
tions at NNLO. This absence is particularly important since NNLO FFNS 
contains terms of the form a^ln (Q'^/mfj). 

At the other extreme we have the ZM-VFNS. Here the terminology 
scheme is misleading. It usually means a different way of arranging the 
complete calculation. In this case there is an overall error of 0{m^/Q'^). 
In my opinion ZM-VFNS is not useful. At high scales we are often in the 
limit where charm and bottom are effectively massless and a GM-VFNS is 
identical to the ZM-VFNS. However, the partons are obtained from fitting 
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Figure 2. Comparison of the CTEQ6M singlet distribution with uncertainties and the 
CTEQ6HQ singlet distribution. 

to data in regions where 0{m?^ IQ"^) corrections are important, and ignoring 
these leads to incorrect partons at aU in the ZM-VFNS. In Fig. 2 we see 
the difference between the GM-VFNS CTEQ6HQ partons^^ ^nd the ZM- 
VFNS CTEQ6 partons with their (conservative) uncertainties. At NNLO 
the partons become discontinuous at the transition points, indeed c{x, Q^) 
at is negative, and at this order we certainly need a GM-VFNS. If 
for some process GM-VFNS coefficient functions are not known, the error 
of 0{rnjj/Q^) from using the GM-VFNS partons is no worse than the 
permanent error from using ZM-VFNS. At worst we can input kinematic 
constraints into coefhcient functions. There are a variety of definitions of 
a GM-VENS"^**, but they generally agree on fundamentals. Each choice 
is superior to ZM-VFNS and to FFNS. However, most are not defined in 
detail up to NNLO^^, and there are some lingering differences. 
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Figure 3. Comparison of the benchmark gluon distributions and dy distributions. 
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For a given theoretical prescription we still have a wide choice of partons. 
It is obvious that some competition is necessary, but not all partons are 
equal - some are, in some sense, incorrect. There are a variety of reasons for 
this - bugs in programs, incorrect theoretical approach (e.g. wrong coupling 
for flavour scheme), approximations to complete theoretical approach, or 
region of applicability, e.g. MRSTOSc partons are only suitable within the 
region of cuts on the data fit. The error is sometimes small, but can be the 
size of the intrinsic uncertainty or greater. If so, such partons should not 
be used. Indeed, NNLO is often still in the approximate stage. 

There is also the issue of the treatment of experimental errors. As an 
exercise for the HERA-LHC^^ workshop, partons were produced from fits 
to HI, ZEUS, and BCDMS^^ structure function data for > 

9GeV^ using ZM-VFNS and a common form of parton inputs at Qq = 
IGeV^ - clearly very conservative. Partons were obtained using the rigorous 
treatment of all systematic errors (labelled Alckhin) and using the simple 
quadratures approach (labelled MRST), both using Ax^ = 1 to define the 
limits of uncertainty. As seen in Fig. 3 there are small differences in the 
central values and similar errors, i.e. the two sets are fairly consistent. Even 
so, the full treatment of systematic errors is presumably better, but perhaps 
it is not so straightforward. Consider the averaged Hl-ZEUS data sets^^, 
where the systematics of one data set can be significantly reduced by fitting 
to the other set. The averaged data set is much more precise with very small 
systematic errors. At the HERA-LHC workshop a comparison was made 
of a fit to both data sets and a fit to the averaged data set"^^. The partons 
resulting can differ by more than the uncertainty in each, and the movement 
of the data relative to the theory was different in each case. Data can move 
relative to theory due to systematic uncertainties, but in reality this may be 
due to failures in theory rather than due to the central values for the data 
being incorrect. The conventional approach to systematic errors assumes 
we fit data to a perfect theoretical model with some unknown parameters, 
whereas in fact we are testing QCD at some order, and it may be slightly 
lacking. It is always best to remember this and try to minimize systematic 
errors. This makes the averaged Hl/ZEUS data set very desirable. It is 
easier to understand and trust dominant statistical errors. 

We should also include as many data as possible in order to determine 
the partons, e.g. we see in Fig. 4 the reduced uncertainty in the ZEUS 
partons when including their own jet data"'^^. The central values can also 
move, but do not do so much in this case. As a more dramatic example 
we consider the HERA-LHC benchmark partons and investigate how these 
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Figure 4. Comparison of the uncertainty on the ZEUS gluon distribution with and 
without the inclusion of their jet cross-section data. 

compare to partons obtained from a global fit (the MRSTOl partons^°), 
where the uncertainty is determined using A^^ = 50. There is an enor- 
mous difference in the central values, sometimes many a, as seen in Fig. 5. 
The uncertainties are similar using Ax^ = 1 compared to Ax^ = 50 with 
approximately twice the data. Moreover, a5(A/|) = 0.1110 ± 0.0015 from 
the benchmark fit compared to asij^l'z) = 0.119 ± 0.002. Something is 
clearly seriously wrong in one of these analyses, and I am very confident 
that it is the benchmark fit. It fails when compared to most data sets 
not included, and not all can be unreliable. Partons should be constrained 
by all possible reliable data. The benchmark fit partons are extreme, but 
some other partons frequently used are similar in terms of the quantity of 
data fit, but many input implicit constraints from elsewhere. Also, for the 
global fit Ax^ = 1 is not reliable. There must be something better than 
Ax^ — 50(100) or the offset method^^, but we are not yet sure what that 
is. The problems are partially due to the strict incompatibility of different 
data sets. Systematic errors are difficult to understand and not usually 
Gaussian in nature. Also our theory is never perfect - it is not simply a 
matter of tying down unknown constants. There are corrections possible at 
low Q^, small x and high x. Indeed, comparing different sets of partons, one 
finds that the gluon is still very uncertain at low x and Q^, even though all 
partons are fit to the same small- a; HERA data. The additional constraint 
from a direct measurement of Fl{x, Q^) would help this situation**^. 
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Figure 5. Comparison of the benchmark gluon distribution and dy distribution with 
the corresponding MRST2001E partons. 

To conclude, there are many types of partons, and although some may 
be ignored, a variety is needed for the full range of applications and to esti- 
mate the uncertainties due to different assumptions (though one should be 
•wary of partons that have only a limited set of constraints from data). We 
need different prescriptions for heavy flavours (though not ZM-VFNS), dif- 
ferent factorization schemes and different orders. As a test of QCD, NNLO 
is preferable, and •we are no'w obtaining reliable NNLO partons. We some- 
times need partons for special occasions, e.g. to investigate the NuTeV 
siri^ 6w anomaly^^. There are many available, -with QED corrections'', 
isospin violation^, s{x,Q^) ^ s{x,Q'^)'^^ etc.. We also need to determine 
whether resummations at small or large a:, higher twist or other theoretical 
corrections are important in some regions. There is much activity in these 
areas, and hopefully it will very soon provide concrete results. 
References 

1. E866 Collaboration: R.S. To'well et al, Phys. Rev. D64 052002 (2001). 

2. NuTeV Collaboration: M. Goncharov et al, Phys. Rev. D64 112006 (2001). 

3. S.J. Brodsky, P. Hoyer, C. Peterson and N. Sakai, Phys. Lett. B93 451 (1980). 

4. A.D. Martin et al, Eur. Phys. ,L C39 155 (2005). 

5. H.L. Lai et al, Phys. Rev. D55 1280 (1997). 

6. H.L. Lai et al. CTEQ Collaboration: Eur. Phys. J. C12 375 (2000). 

7. J Pumplin et al, JHEP 0207 012 (2002). 

8. A.D. Martin et al, Eur. Phys. J. C4 463 (1998). 

9. A.D. Martin et al, Eur. Phys. J. C35 325 (2004). 

10. A.D. Martin et al, Phys. Lett. B604 61 (2004). 

11. S. Alekhin, Phys. Rev. D63 094022 (2001). 

12. S. Alekhin, Phys. Rev. D68 014002 (2003). 

13. M. Cluck, E. Reya and A. Vogt, Eur. Phys. J. C5 461 (1998). 



February 2, 2008 6:27 Proceedings Trim Size: 9in x 6in 



dis06procpartons 



8 

14. W.T. Giele and S. Keller, Phys. Rev. D58 094023 (1998); 
W.T. Giele, S. Keller and D.A. Kosower, hep-ph/0104052. 

15. ZEUS Collaboration: S. Chekanov et al, Phys. Rev. D67 012007 (2003). 

16. ZEUS Collaboration: S. Chekanov et al., Eur. Phys. J. C42 1 (2005). 

17. HI Collaboration: C. Adloff et al., Eur. Phys. J. C21 33 (2001). 

18. M. Botje, Eur. Phys. ,J. C14 285 (2000). 

19. HI Collaboration: C. Adloff et al., Eur. Phys. J. C13 609 (2000); Eur. Phys. 
J. C19 269 (2001); S. Chekanov et al, Eur. Phys. J. C21 443 (2001); Phys. 
Rev. D7G 052001 (2004). 

20. DO Collaboration: B. Abbott et al, Phys. Rev. Lett. 86 1707 (2001); CDF 
Collaboration: T. Affoldcr et al., Phys. Rev. D64 032001 (2001). 

21. A.D. Martin, et al., Phys. Lett. B531 216 (2002). 

22. E.B. Zijlstra and W.L. van Neerven, Phys. Lett. B272 127 (1991); ibid. 
B273 476 (1991); ibid B297 377 (1992); Nucl. Phys. B383 525 (1992). 

23. S. Moch et al., Nucl. Phys. B688 101 (2004); Nucl. Phys. B691 129 (2004). 

24. C. Anastasiou et al., Phys. Rev. Lett. 91, 182002 (2003); Phys. Rev. D69, 
094008 (2004). 

25. N. Kidonakis and J.F. Owens, Phys. Rev. D63 054019 (2001). 

26. R.S. Thornc, Phys. Rev. D73 054019 (2006). 

27. M. Klasen and G. Kramer, Phys. Lett. B386 (1996) 384. 

28. B.W. Harris and J. Smith, Nucl. Phys. B452 109 (1995); Phys. Rev. D57 
2806 (1998). 

29. S. Frixione, P. Nason and B.R. Webber, JHEP 0308 007 (2003). 

30. G. Marchesini and B.R. Webber, Nucl. Phys. B33G 261 (1990). 

31. E. Laenen et al., Nucl. Phys. B392 162 (1993); B.W. Harris and J. Smith, 
Nucl. Phys. B452 109 (1995). 

32. A.D. Martin, W.J. Stirling and R.S. Thorne, Phys. Lett. B636 259 (2006). 

33. S. Kretzer et al, Phys. Rev. D69, 114005 (2004). 

34. M. Aivazis et al, Phys. Rev. D50 3102 (1994); W.K. Tung et al, J. Phys. 
G28 983 (2002); S. Kretzer et al, Phys. Rev. D69 114005 (2004); R.S. Thorne 
and R.G. Roberts, Phys. Lett. B421 303 (1998); Phys. Rev. D57 6871 (1998); 
A. Chuvakin et al, Phys. Rev. D61 (2000) 096004. 

35. M. Dittmar et al, "Parton distributions: Summary report for the HERA - 
LHC workshop," hep-ph/0511119. 

36. NMC Collaboration: M. Arneodo et al, Nucl. Phys. B483 (1997) 3; Nucl. 
Phys. B487 (1997) 3. 

37. BCDMS Collaboration: A.C. Benvonuti et al, Phys. Lett. B223 (1989) 485; 
A.C. Benvonuti et al, Phys. Lett. B236 (1989) 592. 

38. A. Glazov, AIP Conf. Proc. 792, 237 (2005). 

39. A. Cooper-Sarkar and C. Gwenlan, hep-ph/0508304. 

40. A.D. Martin et al, Eur. Phys. J. C28 455 (2003). 

41. R.S. Thorne, the proceedings of "New Trends in HERA Physics 2005" Ring- 
berg Castle, Tegernsee, October 2005, p. 359. hep-ph/0511351. 

42. NuTeV Collaboration: G. P. Zeller et al, Phys. Rev. Lett. 88, 091802 (2002) 
[Erratum-ibid. 90, 239902 (2003)]. 

43. F. Olness et al, Eur. Phys. J. C40, 145 (2005). 



